Resonant scattering of matter wave gap-solitons by optical lattice defects 
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The physical mechanism underlying scattering properties of matter wave gap-solitons by linear 
optical lattice defects is investigated. The occurrence of repeated reflection, transmission and trap- 
ping regions for increasing strengths of an optical lattice defect are shown to be due to impurity 
modes inside the defect potential with chemical potentials and numbers of atoms matching corre- 
sponding quantities of an incoming gap-soliton. For gap-solitons with chemical potentials very close 
to band edges, the number of resonances observed in the scattering coincides with the number of 
bound states which can exist in the defect potential for the given defect strength. The dependence 
of the positions and widths of the transmission resonant on the incoming gap-soliton velocities are 
investigated by means of a defect mode analysis and effective mass theory. The comparisons with 
direct integrations of the Gross-Pitaevskii equation provide a very good agreement confirming the 
correctness of our interpretation. The possibility of multiple resonant transmission through arrays 
of optical lattice defects is also demonstrated. In particular, we show that it is possible to design the 
strength of the defects so to balance the velocity detunings and to allow the resonant transmission 
through a larger number of defects. The possibility of using these results for very precise gap-soliton 
dynamical filters is suggested. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.75.-b 



I. INTRODUCTION 

Bose-Einstein condensates (BECs) in optical lattices 
(OLs) are presently attracting a great deal of interest [l[ 
due both to their flexibility in parameter design and to 
the possibility they offer to observe interesting phenom- 
ena such as superfluid to Mott-insulator transition 0, 
Bloch oscillations Q, Landau-Zener tunneling [^S, gen- 
eration of coherent atomic pulses (atom laser) |6(, atom 
interferometry 0, etc. In this respect, OLs allow to con- 
trol important properties of BEC by means of their peri- 
odic structure, this allowing, for example, the existence 
and stability of localized nonlinear excitations with chem- 
ical potentials inside band-gaps (so called gap-solitons 
(GSs)) even in the presence of repulsive interactions (pos- 
itive scattering lengths). This fact, that would be obvi- 
ously impossible in absence of the OL Q, has been ex- 
perimentally demonstrated in ||. 

Modulations of the OL can be used to accelerate, de- 
celerate or to scatter GSs as well as to control their ve- 
locities [T3|. Uniform accelerations of the OL combined 
with periodic modulations of the scattering length, ei- 
ther in space or in time, were shown to be effective tools 
to induce long-lived Bloch oscillations in the nonlincar 
rcgime (Tlj . as well as band-gap tunneling phenomena 
such as the Landau-Zener tunneling [l2|. Periodic time 
dependent OL accelerations were also used to achieve the 
dynamical localization of nonlinear matter waves (l3l . Il4j 
and the Rabi-oscillations of GS states across a band- 
gaps [l5j which survive on a long time scale in the pres- 
ence of nonlinearity. Moreover, the combination of the 
above phenomena permits the stirring of GSs both in 
the reciprocal and in the direct lattice space as recently 



demonstrated in [lfj. Besides the one-dimensional con- 
texts in which these effects have been investigated, OLs 
also play an important role for stabili zing GSs against 
collapse or decay in higher dimensions [17| ■ 

All what said above refers to the case of perfect OLs, 
e.g. OLs without distortions or defects which compro- 
mise the periodicity structure. In this case the interplay 
between periodicity and nonlinearity is the only source 
for the localization of matter waves in the system. OL 
defects, however, introduce additional bound states in 
the band-gaps (so called impurity modes) (l8| providing 
an alternate source of localization in the system. Oppo- 
site to GSs, impurity modes exist both in absence and in 
presence of nonlinearity and can interfere in the scatter- 
ing process of GSs by OL defects. 

Scattering properties of solitons with single (non peri- 
odic) potential wells have been extensively investigated 
during the past years. For the case of the nonlinear 
Schrodinger equation (NLSE) scattering of solitons by 
extended defects were numerically investigated in [19j 
where the occurrence of a scries of reflection, transmis- 
sion, and trapping regions as a function of the defect 
strength was reported. A similar problem has been re- 
cently investigated for solitons of the Gross-Pitaevskii 
equation (GPE) with rectangular potential wells and for 
attractive interatomic interactions (20j . 

Similar studies were also done for point defects of the 
discrete NLSE equation [2lT - |24l | which corresponds, un- 
der suitable conditions, to a tight-binding model of BEC 
in a deep OL (25j . The resonant transmission, reflection 
and trapping of discrete breathers by point-defects was 
investigated in [lfj . In contrast to these studies, however, 
the scattering of continuous GPE GSs by OL defects have 
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been scarcely investigated. In this context, we mention 
the numerical study performed in [27j where the scat- 
tering properties of GSs were suggested to be useful to 
construct quantum switches and quantum memories. To 
our knowledge, however, the mechanism underlying res- 
onant transmissions of GPE GSs by OL defects and the 
possibility of multiple defects resonant scattering have 
not yet been discussed. 

In the present paper we provide an extensive numerical 
investigation of the scattering properties of GSs by OL 
defects and identify the physical mechanism underlying 
the phenomenon of the resonant transmission. In par- 
ticular, we show that the presence of repeated reflections 
transmission and trapping regions observed for increasing 
strengths of an OL defect is associated to the impurity 
modes inside the defect potential with energies (chem- 
ical potentials) and numbers of atoms matching corre- 
sponding quantities of the incoming GS. As the OL de- 
fect strength is increased, new impurity modes enter from 
the gap edges, moving toward the center of the gap (bot- 
tom of the potential). This implies that for incoming 
GSs with chemical potentials very close to band edges, 
the number of resonances observed in the scattering co- 
incides with the number of bound states which can exist 
in the defect potential for the given defect strength. This 
fact is demonstrated both by studying stationary states 
inside the defect potential and by direct numerical inte- 
grations of the GPE. An excellent agreement between the 
two approaches is found, this confirming the correctness 
of our interpretation. 

The dependence of the resonant transmission on the 
incoming velocity of the GS is also investigated both for 
fundamental GSs in the semi-infinite gap and for GSs in 
the first gap zone. As a result we show that the trans- 
mission resonant peaks become wider as the incoming 
velocity is increased, with very sharp peaks at small ve- 
locities. The multiple resonant transmission through a 
series of (two and three) OL defects is also demonstrated. 
We show in this case that for equally spaced identical OL 
defects the widths of the full transmission resonances in 
general decreases as the number of defects is increased. 
We demonstrate, however, that the resonant transmis- 
sion through a series of defects can be achieved if defects 
are designed so to compensate off-resonance detunings in- 
troduced by velocity changes. This fact gives rise to the 
possibility of using arrays of OL defects as very precise 
filters for matter wave dynamics. 

The paper is organized as follows. In Sec. II we present 
the model equation and discuss the basic properties of 
GSs of the GPE. In the Sec. Ill we consider the interac- 
tion of small amplitude GSs with a localized Gaussian im- 
purity in the OL. The problem is investigated by means 
of direct numerical integrations of the GPE for both at- 
tractive and repulsive interactions as well as for attractive 
and repulsive defects. In section IV we use a stationary 
defect mode analysis to show that the repeated resonant 
transmission, reflection and trapping regions occurs in 
correspondence of resonances with impurity modes inside 



the defect potential and investigate their dependence on 
GS incoming velocities. In section V the multiple res- 
onant transmission across a series of equidistant (equal 
and unequal) OL defects, is investigated. In Sect. VI the 
main results of the paper are shortly resumed. 

II. MODEL EQUATION AND STATIONARY 
LOCALIZED SOLUTIONS 

Let us consider a cigar-shaped BEC described by the 
following normalized one-dimensional GPE 0, [28[ 

= -ipxx + V ext (x)ip + a\ip\ 2 ip, (1) 

where V ex t (x) denotes an external potential of the form: 
V ext (x) = V i(x) + Vd(x), with V i(x) a perfect OL of pe- 
riod L: V i(x) = V i(x + L), and Vd(x) a defect potential 
consisting of a sum of rid single wells potentials localized 
on a distance of several lattice period around the OL sites 
Xi, i = 1, rid- In the following we assume V i and Vd 
to have the form 

V o i(x) = V cos(2x), (2) 

V/ d (z) = £-^-ex P Rz-*,) 2 /(2A.?)] (3) 
i=i v 1 

where without loss of generality the period of the OL is 
taken it. We remark that results of this paper will not 
qualitatively dependent on type of defect (we used Gaus- 
sian defects just for numerical convenience) and similar 
results can be obtained for other shapes of the defects, 
like square well defects, for example). In the following 
we will mainly restrict to the cases of few OL defects: 
rid = 1,2,3, only. We also remark that in Eq.([T]) the 
normalization has been made by measuring the energy 
in units of recoil energy E r = h 2 k 2 /(2m), where k = n/d 
and d is the lattice constant, the space coordinate and 
time in units of d/ir and E r /h, respectively. The dimen- 
sionlcss macroscopic wave function is also normalized as 
J \ip\ 2 dx = 87r./V£;|a s |, where a s is the s-wave scattering 
length. 

It is well known that in the absence of the defect poten- 
tial Vd = 0, Eq.JT]) posses families of exact GS solutions 
with energy located in the band-gaps of the linear eigen- 
value problem 

+ [E a (k)-V ol (x)]tp ak = 0, (4) 

where ip a k (x) are orthonormal set of Bloch functions with 
a denoting the band index and k the crystal-momentum 
inside the first Brillouin zone (BZ): k G [—1,1]. It is 
also known that small-amplitude GSs with chemical po- 
tentials E s very close to band edges are of the form 
ip(x,t) = A(C, r)(p a k(x)e~ lEa ^ t with the envelope func- 
tion A((, t) obeying the following NLSE 
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FIG. 1: Band-gap structure for Vb = —1 and GS bound states 
at points A, B and C near band-gap edges, corresponding 
to chemical potentials Ea,b,c = —0.125; 0.475; 1.452, respec- 
tively. Insets show the wave functions ipA,B,c of the bound 
states. 



where r and £ are slow temporal and spatial variables, 
M e ff = (d 2 E a /dk 2 )~ 1 denotes the soliton effective mass 
and x = a J Wak^dx the effective nonlinearity Q. The 
condition for the existence of such solitons is x^ef / < 
and coincides with the condition for the modulational 
instability of the Bloch wavefunctions at the edges of the 
BZ Examples of small amplitude GSs with chemical 
potential inside the band-gap structure are depicted in 

Fig. m 

In the presence of very diluted OL defects, GSs will 
continue to exist and away from defects they practically 
coincide with GSs of the undistorted OL. An attractive 
(local potential well) or repulsive (local potential barrier) 
OL defect will be seen by the GSs differently, depending 
on the sign of their effective mass. Thus, for example, 
in the case of repulsive interactions and a negative ef- 
fective mass, a GS approaching a repulsive defect will 
see it as a trapping potential (rather than as a potential 
barrier) , thus besides being totally or partially transmit- 
ted/reflected, it can also be trapped at the defect site, a 
fact which would be impossible in absence of the OL. In 
all the numerical simulations presented in this paper we 
have used stationary GSs of Eq s. fB -dSt, exactly deter- 
mined by the shooting method [29[ or by self-consistent 
calculations (30j . and have put them in action by means 
of phase imprinting (e.g. we multiply the state by the 
phase factor e - lcrvx / 2 7 with v being the GS velocity). 
This provides initial condition for the GPE numerical 
time integration of the form: ip(x,0) = i\) s (x)e~ %avx l 2 . 
For possible experimental implementations of our results, 
an alternate method to use to put the GS in action could 
be the acceleration of the OL for a short time interval 
to move the stationary state away from the BZ edges 
(center) so that it can acquire a small Bloch velocity 
Vb = dE„ (k) I dk -C 1 (for details on how this can be 
done see (lo|V 



FIG. 2: Dependence of the transmission (left), trapping (cen- 
ter) and reflection (right) coefficients on the amplitude rj and 
width A of the defect potential Vd- Parameters are v = 0.05, 
E s = -0.125, V = -1. 



III. RESONANT SCATTERING OF GS BY AN 
OL DEFECT 

We consider first the scattering of a GS by a single lo- 
calized defect (n c i = 1 in Eq. ([3])). In order not to perturb 
the soliton initially, the distance between the soliton cen- 
ter and the OL defect is taken much larger than the width 
of the GS (« 8QL). In the following we compute the trap- 
ping, transmission and reflection coefficients defined as 
C = N c {t f )/N ini , T = N T (t f )/N ini and R = 1-(T+C), 

with N C (t f ) = j*i c m f )\Hx, N T (t f ) = mt f )\ 2 dx 

and Ni n i = J_ \ip(t = 0)\ 2 dx denoting the numbers 
of atoms trapped, transmitted and in the initial state, 
respectively. The final time tf depends on the initial ve- 
locity of the GS and in the numerical experiment is deter- 
mined as the time necessary for the coefficients T, R, C 
to become stationary after the scattering process has oc- 
curred. The trapping region [— x c ; x c ] has the size of the 
initial soliton and in all our calculations we fix x c = 30L. 



A. Scattering of semi-infinite GS by an OL defect 

For the scattering of a GS in the semi-infinite gap to 
occur, the existence criterion for a GS near the bottom 
edge of the first band (where the effective mass is posi- 
tive) implies that the nonlinear coefficient must be nega- 
tive sign(er) = — 1. By applying a small initial velocity to 
the GS in the defect direction, depending on the ampli- 
tude r\ and width A of the defect, three possible scenar- 
ios can occur: i) complete reflection, R = 1; ii) complete 
transmission, T = 1; iii) partial trapping, C > 0. The 
regions of the parameter space {77, A} where these dif- 
ferent regimes occur, as obtained from direct numerical 
integration of Eq.([T]), are reported in Fig. [5] The depen- 
dence of T, C, R on the defect strength 77 for two different 
incoming GS velocities and for a fixed value of the defect 
width A = 5 is depicted in Fig. [3] We see that by chang- 
ing the strength of the defect it is possible to achieve 
complete reflections (R w 1) or transmission (T ps 1) as 
well as partial trappings (C > 0). The profiles of the 
initial, reflected and transmitted GS arc depicted in the 
top panels of Fig[3]for defect strengths corresponding to 
points labeled in the middle panel by letters A, B, C. By 
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FIG. 4: Defect modes corresponding to the maxima of the 
trapping curves (C ~ 1) indicated in Fig[3]by letters D, E, 
F and G, respectively. Green lines correspond to the defect 
modes obtained from the dynamical GPE calculations while 
dark lines refer to defect mode profiles obtained from Eq.® 
(see below). Parameters are fixed as in FigO 



FIG. 3: Dependence of coefficients T(rj), R(rj) and C{rj) 
on the OL defect strength for two GS incoming velocities: 
v — 0.01 (middle panel) and v = 0.05 (bottom panel). The 
top panels A, B and C correspond to the profiles of initial, 
reflected and transmitted GS. The wide line in gray shows 
the position of the defect. Other parameters are A = 5, 
E s = -0.125, Vb = -1. 

comparing the middle and the bottom panels of FigJ3] it 
is clear that the sharp peaks at small incoming velocities 
(v fx 0) for which T fx 1 (sec points A, C, E, F, G in the 
middle panel) become wider as the velocity is increased 
while the regions for which R fx 1 are a bit reduced (also 
see FigElfor the case of first band-gap GSs). 

The first four impurity modes corresponding to the 
transmission peaks D,E,F,G in the middle panel of 
Figj3] have been depicted in Fig. [4] Notice the alter- 
nating odd-even symmetry of these modes (modes D and 
F being odd and modes E and G being even with respect 
to the center x — of the defect potential) as usual 
for eigenstates of one-dimensional trapping potentials. 
We remark that the existence of four resonant trans- 
mission peaks (and reflection regions) seen in Fig|3] for 
— 7 < r\ < correlates with the existence of four im- 
purity modes for the given defect strength region (see 
stationary defect mode analysis below). 

B. Scattering of first band-gap GSs by an OL 
defect 

Similar results as those of the previous section can be 
found for GSs inside the first band-gap, with the only 
difference that now there are two possibilities for the ex- 
istence of small-amplitude solitons: i) in the vicinity of 
the top edge of the first band where the effective mass 
is negative and therefore GSs can exist only for repulsive 
interactions a = 1; ii) and in the vicinity of the bottom 



edge of the second band where effective mass is positive 
and GS exist only for attractive interactions a = — 1. 
The small-amplitude GS near the first band-gap edges 
B,C, in FigQ]are shown by the corresponding profiles ipB 
and ipc depicted in the figure. 

As remarked before, the sign of the effective mass de- 
termines the type of the interaction of the GS has with 
the defect potential and for negative GS effective mass 
(repulsive interatomic interactions) the defect will be 
seen as a defect trapping potential (supporting therefore 
bound states) if the defect strength r\ is positive rather 
than negative (as seen for the case of a positive effective 
mass). Except for this, results go in parallel with those 
of the previous section and have been collected in Figs|5J 
[7] for the case of an initial GS close the top edge of the 
first band (negative effective mass). In particular, from 
the bottom panel of Fig. [5] and Fig. [7] we clearly see 
that for a fixed defect width the region of the resonant 
transmission (corresponding to red color) becomes wider 
as the incoming GS velocity is increased, while the full 
reflected regions are achieved mainly for small velocities, 
as one could have expected. We remark that also in this 
case the number of resonant transmission reflection and 
trapping regions are found to correlate with the number 
of impurity modes present in the defect potential for the 
given range of the defect strength (first four modes cor- 
responding to the maxima of the trapping coefficient C 
in FigE] are depicted in FigEl EJ . 



IV. IMPURITY MODE ANALYSIS AND GS 
SCATTERING BY OL DEFECTS 

As remarked before, the presence of a repulsive (at- 
tractive) OL defect in the GPE affects the existing GS 
states for 77 = and introduces additional localized states 
inside the gap in presence of repulsive (attractive) non- 
linearities. 
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FIG. 5: Dependence of the transmission, reflection and trap- 
ping coefficients on the defect amplitude r\ for A = 5 and 
for two different incoming velocities of the GS (E a = 0.475): 
v — 0.02 (middle panel) and v — 0.1 (bottom panel). Top 
panels A-C show the initial, reflected and transmitted GS 
profiles, respectively. 



Numerical calculations show that the band structure is 
only slightly affected by an OL impurity, the main effect 
being the introduction of bound states spatially local- 
ized at the impurity sites and with chemical potentials 
inside the band-gap. Note that repulsive (attractive) OL 
defects in the presence of an attractive (repulsive) nonlin- 
earity cannot introduce additional bound states because 
the corresponding impurity potentials correspond to bar- 
riers rather than potential wells, due to the positive (neg- 
ative) effective mass. Recalling that the opposite signs of 
nonlincarity and effective mass is a necessary condition 
for the GS existence, one has that the effective impu- 
rity potential acts as a trapping potential when r\ and a 
have equal signs. Away from the OL impurity localized 
states are practically the same as for rj = case. At 
the defect site, however, GSs levels get slightly shifted 
by the impurity potential and additional impurity modes 
enter the gap. A GS moving through the impurity will 
have, in general, a mismatch in energy and in number 
of atoms with the impurity modes, this giving a partial 
reflection/transmission of the incoming matter wave. 

A total transparency of the OL defect is expected 
for incoming GS energies and number of atoms exactly 
matching those of an impurity modes inside the defect 
potential (notice that the energy of impurity modes de- 
pend on rj and on the number of atoms). By increas- 
ing the strength of the impurity, the depth of the defect 
potential increases and more impurity modes enter the 
gap. As |?7| is increased, the energies of these modes en- 
ter the gap from the top (bottom) of a band for impurity 
strength and the nonlinearity both positive (negative). 
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FIG. 6: First four defect modes (panels E-H) in correspon- 
dence of maxima of the C(ri) curve in Fig. [5] Parameters are 
the same as in Fig. [5] 



This implies that the transparency (complete transmis- 
sion) of the impurity occurs in correspondence of each 
impurity mode entering the gap and matching the en- 
ergy of the incoming GS given by 

1 



E 



2 

.(0) 



eff V B 



(0) 



where the first term Eg represents the energy of the 
stationary GS state at k = fco (eg. with ve(fco) = 0), 
while the second one is the contribution due to the kinetic 
energy (here v B (k) = dE a (k)/dk is the Bloch velocity 
and M e ff = (d 2 E a j dk 2 )~ 1 \k=k a the effective mass). For 
small velocities the energy (j6)) practically coincides with 

Es°^ but in general the kinetic term should be accounted 
in the matching in energy with the impurity levels (see 
below). Notice that Eq. ^ is only valid near stationary 
points fco (bottoms or tops of a band) where vs{k) « 
(fc — ko)/M e ff and E can be written as 



EaE p + HLZ^l. (7) 

In the range of initial velocities vb = [0,0.2] we have 
considered, the energy curves E\{k) in vicinity of fco = 
and fco = 1 are very well approximated by Eq. ([7]) with 
£i(0) = -0.12177, £>i(l) = 0.47065, and M eff « 0.565 
and M e ff ss —0.167 for bottom and top edges of the 
band, respectively. By knowing vb (e.g. fc) and M e ff one 
can compute the energy shift due to the nonzero velocity 
to be accounted in the matching between the GS and the 
impurity levels (see lower panels of Figs. |9]). Notice 
that the kinetic energy has the sign of the effective mass 
so that E is pushed forward the corresponding band edge 
for finite vb, meaning that inside the impurity potential 
the GS matching condition with an impurity mode can 
be achieved for a lower values of \rj\. From this we expect 
the resonance transmission peaks to be shifted away from 
the v = resonance toward lower values of \rj\ as vb is 
increased. 

To confirm these prediction with GPE calculations we 
have solved the stationary problem 

u xx + [E -Voi(x) -V d {x)]u-au 3 = (8) 

searching for bound states u(x) localized at the OL defect 
using both shooting method and self-consistent calcula- 
tions. Using these approaches we found that the values 
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FIG. 7: Dependence of C(r],v) (top panel) and T(rj, v) 
R(?7, v) (bottom panel) on r\. Parameters are fixed as: E s = 
0.475, A = 5, V = -1, a = 1. 



of 77 for which the Gaussian defect becomes transparent 
to the GS dynamics, correspond to the values for which a 
localized mode inside the defect potential has energy and 
number of atoms matching the corresponding quantities 
of the incoming GS. Since the discussion is very similar 
(a part implications due to signs of the effective masses) 
for GS of the semi-infinite and for the ones inside finite- 
gaps, we refer to the case of an initial GS inside the semi- 
infinite gap, with a positive effective mass with chemical 
potential close to the bottom of the first band. In the top 
panel of Fig|5]we have shown the dependence of the num- 
ber of atoms N in a defect mode of a given symmetry on 
77, for an energy E = E s = —0.125 corresponding to a GS 
with zero incoming velocity (point A of Figs, [l] [3J . The 
horizontal dotted line refers to the numbers of atoms in 
the initial GS. In the bottom panel of Fig. [5]we show the 
energy mismatch \E — E s \ between the energy of a defect 
modes E and the one of an incoming GS given by Eq. 
[7] for different incoming velocities. Notice that the inter- 
section points A,D,E,F,G of the dotted line N = Na 
with the curves N(rj) are in coincidence with the zeros of 
the function \E — E s \ displayed in the bottom panel for 
v ~ and correspond to the maxima of the transmission 
coefficient T(rf) in Figj3] (see curves v ps in the middle 
panel). Fig. [5] (see also |H] for repulsive case) also explains 
the decay of the reflection coefficient and the rising of 
the trapping coefficient coinciding up to the maxima ob- 
served (just before the resonance) for increasing values 
of \rj\ away from the resonance points A, D, E, F, G. To 
the right of these points in FigJSl the number of particles 
in the defect mode is higher then the number of particles 
in the incoming soliton N(rj) > Na- In this case the GS 
cannot "use" a defect mode to pass the OL impurity and 
will be fully reflected. On the left of the resonance points, 
the number of atoms in the GS is higher than the one in 
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FIG. 8: Top panel. Number of atoms in defect modes vs 
the amplitude of the defect 77. The (black,solid)/(red, dashed) 
lines correspond to the even/odd symmetry of the defect 
modes. Other parameters are E s — —0.125, A = 5, Vo = — 1, 
a = — 1. Bottom panel. Mismatch between defect and GS en- 
ergies, \E — E s \, vs 77 for different incoming velocities: v ~ 
(gray, •), v — 0.05 (red, ■) and v = 0.1 (blue, +). 



the defect mode so that the GS can be captured by the 
OL impurity by releasing the excess number of atoms into 
the reflection and transmission channels. It is clear that 
the peak of the capture coefficient occurs just before the 
resonant transmission (e.g. for < Na — N(rj) <C 1) of 
the GS through the defect (achieved when the condition 
N(rf) = Na is exactly fulfilled). Notice from FigJUthat 
at the resonant transmission, the profiles of the station- 
ary impurity modes obtained by solving Eq.® exactly 
coincide with the solution of the GPE during the passage 
trough the defect. In Fig[TU]we have also depicted, simi- 
larly to FigJ^J the level curves TV = Na as a function of 
the amplitude 77 and the width A of the defect. From 
the bottom panel of Fig. [5] it is also clear that for v 7^ 
the resonances shift in the direction of lower values of 
1 77 1 (as expected from our analysis). This well correlates 
with the GPE calculations reported in the bottom panel 
of FigJ3] Similar results are obtained for a GS near the 
bottom of the first band gap (point B in Fig. [TJ for the 
case of repulsive interactions. This is shown in Fig|9l 
from which we see that the shift of the resonance due to 
the finite Bloch incoming velocity is always in the direc- 
tion of lower values of rj and are in good agreement with 
the GPE numerical results reported in the bottom panel 
of Fig. [SJ In particular, for the velocity vb =0.1 and 
the resonance near r\ = 4 we obtain from Fig. |H]that the 
resonance peak is at 77 rj 4.45 while from the GPE result 
in Fig. [5]we obtain the value 77 ps 4.3. It is also worth to 
note from the bottom panels of Figs. [HI [HI that for a fixed 
energy mismatch the widths of the curves increase as the 
incoming velocity is increased, a fact which correlates 
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FIG. 9: The same as in Fig[8]but for E a = 0.475 and a = 1. 
The incoming soliton velocities are u « (gray, •), v — 0.1 
(red, ■) and v = 0.2 (blue, +). 



FIG. 11: Transmission, reflection and trapping diagram of 
GSs with velocities v = 0.01 (a) and v = 0.05 (b) in presence 
of two Gaussian defects of the OL placed at positions x — 
and x = 100-7T. Other parameters are the same as in Fig[3] 



with the broadening of the transmission peaks observed 
in the GPE calculations (compare bottom panels in Fig. 
[3] and in Fig [SI respectively). 

From this we conclude that the above impurity mode 
analysis fully confirms the interpretation of the trans- 
mission peaks as resonances between incoming GSs and 
impurity modes. 
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FIG. 10: Level curves N = Na for even (red, dashed) and 
odd (black, solid) defect modes vs (77, A) for the case of semi- 
infinite gap (a) and bottom of the first gap (b). Other param- 
eters are fixed as E s = -0.125(0.475), Vo = -1, a = -1(1). 



V. RESONANT TRANSMISSION THROUGH 
MULTIPLE OL DEFECTS 

An interesting question to ask is whether the soliton 
resonant transmission could survive multiple impurity 
scatterings. In Fig. [TT] we show the transmission, re- 
flection and trapping regions of an attractive GS of the 



GPE with two identical Gaussian OL impurities, placed 
at x = and x = 1007T, respectively. The TCR curves de- 
pend also on the distance between the two impurities and 
this could be varied so to find optimal values for double 
resonant transmission to occur. The phenomenon, how- 
ever, is more sensitive to variations of the incoming GS 
velocity as one can see by comparing Fig llll with the case 
of single impurity in Fig. [3] Notice that while the trans- 




FIG. 12: Contour plot of the GS space-time dynamics for an 
OL with two equal Gaussian defects of strengths: 77 = —0.43 
(top left), 77 = -0.57 (top right), 77 = -0.58 (bottom left), 
77 = —0.68 (bottom right). The initial velocity of the GSs is 
v = 0.05 and parameters are the same as for Fig. 1111 (bottom 
panel) . 

mission resonance peaks at small velocities are practically 
unaffected by the presence of the second impurity, they 
become more narrow at larger velocities for the resonant 
scattering on two impurities. The shrinking of the reso- 
nances at larger velocities can be understood by the fact 
that for higher incoming velocities the velocity of the GS 
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after the first impurity is slightly reduced and the vari- 
ation introduces a detuning from resonance in the scat- 
tering with the second impurity which in turn reduces 
the double resonance width. From this one can expect 
that in the presence of more lattice defects the multiple 
resonance transmission peaks become very narrow and 
only solitons with very precise initial matching velocities 
will be able to pass. In general the GS may be able to 
pass only a finite number of impurities before remaining 
trapped at one impurity or becoming scattered back and 
forth between them. This is shown in Fig |12l where con- 
tour plots of the GS space-time dynamics in presence of 
two identical Gaussian defects are reported for different 
strengths 77 and for the same parameters as for the bot- 
tom panel of FigJTTJ In the left top and bottom panels 
we see the occurrence of total reflection and transmission 
for a resonant value of rj while in the corresponding right 
panels we show the back and forth dynamics of the soliton 
between two impurities (top panel) and the trapping by 
the second impurity (top panel) for a non resonant values 
of 77. Notice the small change of the soliton velocity after 
the passage through the first defect. To overcome the de- 
tuning from resonance induced by the change of velocity 
one could design the strength of the second impurity so 
to match the value of the intermediate velocity and still 
achieving T m 1. In this way one can realize a double 
(or multiple) filtering of the soliton motion so that only 
GSs having very precise initial velocity can overcome the 
defect series. This possibility is illustrated in Fig. [13] 
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FIG. 13: Panel (a). Contour plot of the GS space-time dy- 
namics in an OL with three equal Gaussian defects at posi- 
tions Xi — 0, X2 = 60-7T, xs = 120-7T of strengths 771 = 772 = 
773 = —0.67. Panel (b). Case of three different Gaussian de- 
fects of strengths 771 = —0.67, 772 = 771 — 0.02, 773 = 771 — 0.04 
(b). Other parameters are fixed as in Fig[3] 

where the resonant transmission of a GS though a series 
of three Gaussian defects of the OL, with defect strengths 
designed so to compensate the detunings introduced by 
the velocity changes, is shown. Notice from the left panel 
of Fig. [T3] that in the case of equal defects with the ini- 
tial velocity matching the resonance transmission peak 
of the first defect, the GS cannot be transmitted through 
all the three defects but it is stopped at the second de- 
fect. The dynamics of the soliton in presence of multiple 



defects may be quite complicated and is beyond the aim 
of this work (an investigation in the parameter space of 
the possible scenarios for the GS time evolution will be 
discussed elsewhere). 



VI. CONCLUSION 

In this paper we have investigated the scattering prop- 
erties of matter wave gap-solitons with optical lattice de- 
fects in the framework of the mean-field Gross-Pitaevskii 
equation. We have shown that the occurrence of repeated 
reflection, transmission and trapping regions are in corre- 
spondence of the defects strengths for which the number 
of atoms and energies of additional bound states created 
by the optical lattice defect, match the ones of the incom- 
ing gap-soliton. This has been demonstrated by a study 
of the stationary defect modes energies (chemical poten- 
tials) and number of atoms as a function of the defect 
strength. A very good agreement between the predicted 
values of the resonant transmission peaks by means of 
impurity modes and the ones found by direct time inte- 
grations of the Gross-Pitaevskii equation, is found. 

The behavior of the reflection and trapping curves have 
also been explained by impurity mode analysis. These in- 
vestigations have been performed both for attractive and 
repulsive interactions and for localized states both in the 
semi-infinite gap (attractive case) and in the first gap 
zone (attractive and repulsive cases). The dependence of 
the resonant transmission on the initial gap-soliton ve- 
locity has been also investigated. We have shown that 
the positions of resonant transmission peaks shift toward 
lower values |»7| as the initial gap-soliton velocity is in- 
creased, while the widths of resonances shrinks to zero at 
very small gap-soliton velocities. The possibility of mul- 
tiple resonant transmission through an arrays of defects 
was also demonstrated. In particular, we have shown 
that for an optical lattice with two equally spaced iden- 
tical Gaussian defects the widths of the full transmission 
resonances for larger incoming velocities, decreases as the 
number of defects is increased. The sharpening of the 
transmission peaks has been explained in terms of the 
detuning from resonance introduced by the small veloc- 
ity change after the passage of an optical lattice defect. 
Finally, the resonant transmission of a gap-soliton though 
a series of Gaussian defects with unequal strengths de- 
signed so to compensate detunings introduced by the ve- 
locity changes, was demonstrated. These results give the 
possibility to construct very precise filters for the matter 
wave gap-soliton dynamics by means of properly designed 
arrays of OL defects. 
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